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Increased adipocyte size and number are associated with many of the adverse effects observed in metabolic disease 
states. While methods to quantify such changes in the adipocyte are of scientific and clinical interest, manual methods to 
determine adipocyte size are both laboriousand intractable to large scale investigations. Moreover, existing computational 
methods are not fully automated. We, therefore, developed a novel automatic method to provide accurate measurements 
of the cross-sectional area of adipocytes in histological sections, allowing rapid high-throughput quantification of fat cell 
size and number. Photomicrographs of H&E-stained paraffin sections of murine gonadal adipose were transformed using 
standard image processing/analysis algorithms to reduce background and enhance edge-detection. This allowed the 
isolation of individual adipocytes from which their area could be calculated. Performance was compared with manual 
measurements made from the same images, in which adipocyte area was calculated from estimates of the major and 
minor axes of individual adipocytes. Both methods identified an increase in mean adipocyte size in a murine model 
of obesity, with good concordance, although the calculation used to identify cell area from manual measurements 
was found to consistently over-estimate cell size. Here we report an accurate method to determine adipocyte area in 
histological sections that provides a considerable time saving over manual methods. 



Introduction 

Unbiased approaches to the quantification of biological features 
in histological sections will facilitate high-throughput histomor- 
phometry. A number of methods to assess adipocyte size and 
distribution have been described,' "^ but these do not provide the 
flexibility necessary for large-scale studies. Adipocyte size reflects 
adipocyte physiology, increased size being associated with com- 
promised metabolism and the detrimental effects of obesity and 
type 2 diabetes.^'^ Adipose mass can expand by two mechanisms: 
hyperplasia (an increase in cell number) and hypertrophy (an 
increase in cell size).' Automated methods would ideally deter- 
mine both phenomena, allowing, for example, the effects of 
genetic manipulation or therapeutic intervention to be evalu- 
ated. Traditionally, individual adipocytes are measured micro- 
scopically following coUagenase digestion of fat pads; however, 
this does not lend itself to retrospective or large-scale analysis 
and may lead to the loss of larger adipocytes.^ Furthermore, this 
method utilizes equipment not available to all laboratories and 
the osmium tetroxide protocol may cause cell swelling result- 
ing in an overestimation of adipocyte size.^ Manual techniques 
to measure the cross-sectional area of adipocytes in histologi- 
cal sections were subsequently described, but these are labori- 
ous and demand discrimination, counting, and measurement 
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of many hundreds of irregular events.'*'' This limits the num- 
ber of fields that may be accurately counted in a timely fash- 
ion. More recently, manual quantification has been superseded 
by computational methods. ''^'"^ Flowever, some of those meth- 
ods described to date are not fully automated as significant user 
input is required; thus the influence of human error cannot be 
entirely ruled out and further errors may be introduced in some 
of these methods as the formula used to calculate cellular area 
provides only an approximation. The technical issues facing the 
implementation of a fully- automated method are manifold as 
variation in staining intensities, tissue quality (particularly the 
fragmentation of the fragile adipocyte cell wall) and photomi- 
crographic capture must be addressed. Moreover, an automated 
method must deploy a consistent set of rules in discriminating 
cells from artifacts, in order to minimize errors over a large num- 
ber of images. 

Image analysis is a multistep process. Original images must 
be pre-processed to improve signal-to-noise ratios and enhance 
feature segmentation and edge detection. Post-processing algo- 
rithms to deal with, for example, broken cells walls are essential if 
accurate measurements are to be made. Finally, the data must be 
output in a meaningful way. The method described herein uses 
existing principles in image analysis to provide a rapid automated 
determination of adipocyte size from H&E-stained histological 
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Figure 1. Manual adipocyte measurement. Original H&E image of gonadal fat tissue (A). Manual annotation showing minor and major adipocyte axes 
(B). Diagram contrasting estimated elliptical area (black) to actual cell area (red) (C). 



sections. One major advantage of this technique is that it returns 
accurate cell area information as compared with manual meth- 
ods. By comparison of our algorithm to a manual method, we 
were able to confirm its accuracy on adipocyte quantification in a 
murine model of obesity. 

Results 

Whole-slide images were captured using a digital pathology plat- 
form (Aperio ScansScope CS). Using at least 16 representative 
fields of view per whole-slide image (a typical field of 160000 
|jLm~ is shown in Fig. 1 A) manual measurements of adipocyte size 
were calculated. The perpendicular maximum (a) and minimum 
axes (b) were estimated for each adipocyte (Fig. IB) by annotat- 
ing images using Aperio ImageScope software. Adipocyte size 
(|jLm-) was approximated as the area of the ellipse whose major 
and minor axes are a and b, respectively, and calculated using the 
following equation (Fig. IC): 
Ellipse area = T;ab 

These images were then subjected to analysis with a fully-auto- 
mated method. The saturation channel of the hue, saturation, 
value (HSV) color space clearly revealed adipocyte boundaries 
(Fig. 2A and B) and adipocytes were counted after extending 
damaged boundaries and excluding incomplete cells located at 
the edge of each image (Fig. 2C). The absolute area in pixels of 
each object was calculated and converted to |xm- (Fig. 2D). 



Both methods identified a significant {P < 0.01) increase in 
mean adipocyte size in tissue isolated from offspring of high- 
fat fed dams as compared with the offspring of chow-fed con- 
trols (Fig. 3). However, we found that the automated method 
consistently returned smaller values for cell size, the difference 
being exacerbated in larger cells. Specifically, the manual method 
returned an average fat cell size of 1004 |JLm^ in lean animals 
and 1733 |JLm^ in the high-fat group, compared with 901 |xm^ 
and 1350 jJtm^ respectively by automated counting. This may 
have been as a consequence of over-estimates of area achieved 
with manual methods that rely on the calculation of an ellipse 
to approximate cell area. Specifically, ellipses are calculated from 
the longest axis values (which are themselves based in human 
judgment), as demonstrated in Figure IC, and thus the area value 
obtained will consistently be a marginal over-estimate for a cell 
of irregular shape. In turn, this may explain the slight decrease in 
significance observed when diet-associated changes in fat size were 
assessed by the automated method compared with the manual 
approach. Similar changes in size distribution were observed in 
the offspring of high-fat fed groups with both methods (Fig. 4A 
and B). When the graphs of the two methods are superimposed, 
a slight leftward shift in the automated value is seen (Fig. 4C and 
D), consistent with more accurate size measurement as described 
above. The only point where a significant difference {P < 0.0001) 
was observed between automated and manual quantification 
(using a two-way ANOVA with a Bonferroni post-hoc test) was 
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Figure 2. Automated adipocyte size measurement. (A) Original H&E-stained image. (B) Saturation channel isolated from HSV color space. (C)Binary 
image following intensity and morphological transformations. (D) Automated annotation. 



in the smallest size category (<500 |xm^) in the offspring of the 
high-fat fed animals (Fig. 4D) . 

Discussion 

We have developed a novel automatic method to provide accu- 
rate measurements of the cross-sectional area of adipocytes in 
histological sections. Performance of our algorithm was consis- 
tent with routinely utilized manual methods'" " and, moreover, 
we were able to establish increased accuracy in unsupervised 
area measurements. More specifically, we have determined 
that the manual method of calculating the cross-sectional area 



of adipocytes uses assumptions that lead to an over-estimation 
of cell size. The calculation assumes that adipocytes possess a 
regular shape, and that manual axis measurements represent 
the true maximal depth and width for any cell. However, we 
have demonstrated that this leads to an over-estimation of cell 
size by calculating the area of pixels within the cell boundary. 
Although the overall concordance between the two methods was 
high, we observed a statistical difference in the lowest cell size 
category (<500 |JLm^). We suggest that this is as a consequence of 
smaller adipocytes in the high-fat animals being more irregular 
in shape compared with the chow fed group, which is exacerbated 
by more arbitrary exclusion criteria in manual compared with 
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computational counting, leading to a 
slight over-representation of adipocytes 
in the smallest category. 

One problem facing both super- 
vised and unsupervised methods is in 
determining if very large cells are truly 
individual objects, or arise from a loss 
or perturbation in separating plasma 
membranes. We purposefully avoided 
using a watershed-based method that 
splits large cells according to arbitrary 
criteria, such as that used by the freely 
available Image] Adipocytes Tool.'' We 
find that such methods over-estimate 
cells in the smallest size category. We 
are confident that the size distribution 
histograms presented in this report do 
not indicate an overrepresentation of 
cells in the largest category (which will 
contain any conflated groups of cells) , 

particularly as we took steps to minimize processing artifacts 
during the preparation of tissue sections. Algorithmic methods 
to overcome problems in the analysis of histological sections with 
extensive fracturing of cell membranes will necessarily yield inac- 
curate results. 

Overall we suggest that unlike manual counting and exist- 
ing computational approaches, the method described herein can 
provide whoUy-unsupervised, fast, high-throughput batch quan- 
tification of fat morphology in different states. 

Material and Methods 

Animals. Breeding, housing, and procedures were conducted 
in accordance with the UK Government Animals Scientific 
Procedures Act 1986 and approved by The University of 
Buckingham Ethical Review Board. Female offspring used in 
this study were generated in house from breedings of C57B16/J 
dams (Charles River) fed either standard laboratory chow con- 
taining 10% fat by energy (BKOOIE, B&K Universal Ltd) or a 
high-fat diet providing 42% calories by energy (829100, SDS 
Diets) throughout gestation. All mice were sacrificed at 6 weeks 
of age using a schedule 1 method and white gonadal adipose 
depots dissected for histological analysis. 

Histology and image capture. Two tissues per animal (chow, 
« = 4 animals; high-fat, « = 5 animals) were fixed in 10% neu- 
tral buffered formalin for 6 h and dehydrated as standard 
before embedding in paraffin wax. Sections (4 |jLm) were cut 
and mounted on positively-charged glass slides and hematoxy- 
lin and eosin (H&E) staining was performed as standard. All 
tissues described in this study were processed simultaneously 
to minimize artifacts. Bright-field whole-slide images were cap- 
tured at 40 X optical magnification with an Aperio ScanScope CS 
instrument. 

Automated size measurements. Image analysis algorithms 
were applied to the same images used for manual calculations, 
and implemented in MATLAB (R2011b, MathWorks). Images 




Figure 3. Comparison of mean adipocyte size derived from manual and automatic methods. Mean 
adipocyte size ± SD was calculated using at least 16 fields per whole-slide image, from at least 
4 animals per group, chow (CH) of high-fat (HF) fed, using both manual ellipse measurements (A), and 
automatic pixel counting (B). 



were transformed from red, green blue (RGB) color space into 
HSV color space. Automatic background subtraction was per- 
formed on the saturation channel only with a median filter of 
71 X 71 pixels. Feature enhancement (to improve subsequent 
edge detection) was achieved with a Piecewise Linear Grey Level 
Transform (also known as gray level slicing) to enhance a speci- 
fied set of gray level values.'^ The resulting image was binarized 
to allow the application of standard morphological operations. 
First an open operator using a disc structuring element was used 
to separate connected cells by breaking connecting plasma mem- 
branes. To maintain the integrity of the cell membranes, a dila- 
tion operator with another disk structuring element was used, 
which helps to re-enforce boundary pixels eroded by the opening 
operator. Finally, the MATLAB function imclearborder (image) 
was used to clear borders and so exclude incomplete cells located 
at the edge of the image to generate the final transformed image. 
Adipocytes were then counted, and the absolute pixel area of each 
object was calculated and converted to |JLm-. Just as a human 
operator will exclude structures that are too small to accurately 
classify as cells, our algorithm was designed to exclude objects 
<240 jJtm^, which may include cross-sections of fibroblastic 
preadipocytes,'^ as artifacts. A fully functional copy version of 
our software is available for download (with example images) at 
http://webspace.buckingham.ac.uk/klanglands/. 

Statistical analysis. The number and size of adipocytes per 
each field of view were calculated and then summed to give val- 
ues per animal comparing both manual and automatic methods. 
From these values either the mean size was calculated, or size dis- 
tributions were plotted using 10 bins: 1 = 240—499 |JLm~; 2 = 
500-999 |JLm^; 3 = 1000-1499 |xm^; 4 = 1500-1999 fJim^; 5 = 
2000-2499 |JLm^; 6 = 2500-2999 jJim-; 7 = 3000-3499 |xm^; 8 = 
3500-3999 |xm^; 9 = 4000-4499 fJim^; 10 > 4500 jjuml Data 
derived from at least four different animals is presented per chow 
or high-fat maternal group. Pair-wise comparisons were per- 
formed using the Student t test, and distributions by 2-way (diet; 
size category) ANOVA in Prism 5 (GraphPad Software Inc.). 
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Figure 4. Comparison of adipocyte size distributions created using manual and automatic methods. Comparison between chow (CH) and high-fat 
(HF) fed adipocyte size is shown for manual (A) and automated methods (B). Cell size distributions in depots from CH animals were very similar when 
either manual or automated methods were used (C). Similar results were also seen for measurements from HF-diet fed animals, although an underesti- 
mation of the smallest groups (<500 jj-m^) was observed with manual counting, as was a slight rightward shift in the manual curve, consistent with size 
over-estimation (D). 
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